// aoz using raw fit parameters for x_F1 and x_ETF

void fitAoZ1(){
  const double p[] = {
    2.04378e+00, 2.09528e+00, 2.13842e+00, // 20-Ca40,41,42,43
    2.04793e+00, 2.09290e+00, 2.13863e+00, 2.18882e+00, // 21-Sc42,43,44,45
    2.04904e+00, 2.09115e+00, 2.13406e+00, 2.18048e+00, // 22-Ti44,45,46,47
    2.04785e+00, 2.08999e+00, 2.13084e+00, 2.17277e+00, // 23-V46,47,48,49
    2.05002e+00, 2.08730e+00, 2.12755e+00, 2.16638e+00, // 24-Cr48,49,50,51
    2.05013e+00, 2.08657e+00, 2.12449e+00, 2.16346e+00, // 25-Mn50,51,52,53
    2.08766e+00, 2.11811e+00, 2.15181e+00, // 28-Ni57,58,59
    2.09126e+00, 2.11498e+00, 2.14697e+00, // 30-Zn61,62,63
    2.11527e+00, 2.14298e+00, // 31-Ga64,65
    2.11404e+00, 2.14035e+00, // 33-As68,69
    // 2.12054e+00, 2.13968e+00, // 35-Br72,73
    // 2.12507e+00, 2.14192e+00, 2.23959e+00, // 36-Kr74,75,78
    // 2.14022e+00 // 37-Rb76
  };
  static const int n = sizeof(p)/sizeof(double);
  double z[n] = {
    40/20., 41/20., 42/20., 43/20., // 20-Ca40,41,42,43
    42/21., 43/21., 44/21., 45/21., // 21-Sc42,43,44,45
    44/22., 45/22., 46/22., 47/22., // 22-Ti44,45,46,47
    46/23., 47/23., 48/23., 49/23., // 23-V46,47,48,49
    48/24., 49/24., 50/24., 51/24., // 24-Cr48,49,50,51
    50/25., 51/25., 52/25., 53/25., // 25-Mn50,51,52,53
    57/28., 58/28., 59/28., // 28-Ni57,58,59
    61/30., 62/30., 63/30., // 30-Zn61,62,63
    64/31., 65/31., // 31-Ga64,65
    68/33., 69/33., // 33-As68,69
    // 72/35., 73/35., // 35-Br72,73
    // 74/36., 75/36., 78/36., // 36-Kr74,75,78
    // 76/37. // 37-Rb76
  };

  TGraph *g = new TGraph(n, p, z);
  g->SetMarkerColor(2);
  g->SetMarkerStyle(20);
  TCanvas *c = new TCanvas("c", "calibration of Z", 800, 600);
  g->Draw();
  g->Fit("pol1");
} // end function fitZ
